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ABSTRACT 

In this Letter, we suggest that the relativistic protons powering the out- 
flows emanating from radio-loud systems containing black holes are accelerated 
at standing, centrifugally-supported shocks in hot, advection-dominated accre- 
tion disks. Such disks are ideal sites for first-order Fermi acceleration at shocks 
because the gas is tenuous, and consequently the mean free path for particle- 
particle collisions generally exceeds the thickness of the disk. The accelerated 
particles are therefore able to avoid thermalization, and as a result a small frac- 
tion of them achieve very high energies and escape from the disk. In our approach 
the hydrodynamics and the particle acceleration are coupled and the solutions 
are obtained self-consistently based on a rigorous mathematical treatment. The 
theoretical analysis of the particle transport parallels the early studies of cosmic- 
ray acceleration in supernova shock waves. We find that particle acceleration 
in the vicinity of the shock can extract enough energy to power a relativistic 
jet. Using physical parameters appropriate for M87 and Sgr A*, we confirm that 
the jet kinetic luminosities predicted by the theory agree with the observational 
estimates. 

Subject headings: accretion, accretion disks - - hydrodynamics - - black hole 
physics — galaxies: jets 

1. INTRODUCTION 

Radio-loud active galactic nuclei (AGNs) are thought to contain supermassive cen- 
tral black holes surrounded by hot, two-temperature, advection-dominated accretion flows 
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(ADAFs) with significantly sub-Eddington accretion rates. In these disks, the ion tempera- 
ture Tj ~ 10 12 K greatly exceeds the electron temperature T e ~ 10 10 K (e.g., Narayan, Kato, 
& Honma 1997; Becker & Le 2003). The observed correlation between high radio luminosi- 
ties and the presence of outflows suggests that the hot ADAF disks are efficient sites for 
the acceleration of the relativistic particles powering the jets. This motivated Blandford & 
Begelman (1999) to consider the possibility of strong outflows of matter and energy from such 
systems, under the assumption of self-similarity and Newtonian gravity. The same scenario 
was extended to the general relativistic case by Becker, Subramanian, & Kazanas (2001). 
However, a single, global, self-consistent model for the structure of the disk/jet that includes 
a specific microphysical acceleration mechanism has not yet appeared in the literature. In 
this Letter, we explore the consequences of the presence of a shock in an ADAF disk for the 
acceleration of the observed nonthermal jet particles. 

It is well known that for certain sets of outer boundary conditions describing the specific 
angular momentum and the specific energy of the gas supplied to the accretion flow, both 
shocked and shock-free solutions for the global disk structure are possible, although there is 
an upper limit on the amount of viscosity beyond which no shocks can form. However, it is 
clear that shocks can exist in both dissipative and inviscid flows (e.g., Chakrabarti & Das 
2004). When both shocked and shock- free solutions are possible for the same set of upstream 
conditions, one must incorporate additional physical principles in order to determine which 
flow solution actually occurs. Since in general the shock solution possesses a higher entropy 
content than the shock-free solution, we argue based on the second law of thermodynamics 
that shocks are preferred (see Becker & Kazanas 2001). Particle acceleration has not been 
studied in shocked disks in either the viscid or the inviscid cases, and therefore we concentrate 
here on the simpler situation of inviscid (adiabatic) flow. In this Letter we focus on disks 
containing isothermal shocks because the jump in the energy flux at such a shock can be 
identified with the energy carried away from the disk by the relativistic particles in the 
outflow. 



2. DYNAMICAL MODEL 

The dynamical structure of the accretion flows considered here follows the model studied 
by Chakrabarti (1989a) and Abramowicz & Chakrabarti (1990), which comprises an inviscid 
disk that either contains an isothermal shock, or else is shock-free. The model is based on the 
standard one-dimensional conservation equations for advect ion-dominated accretion flows, 
including utilization of the pseudo-Newtonian gravitational potential to simulate the effects 
of general relativity (Paczyhski & Wiita 1980). In general, the flow is assumed to be adiabatic 
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and in vertical hydrostatic equilibrium, although the entropy can jump at the shock. Under 
the inviscid flow assumption, the disk/shock model depends on two parameters, namely the 
energy transport rate per unit mass, e, and the specific angular momentum, £. The value of 
e is constant in ADAF type flows since there are no radiative losses, although it will jump at 
the location of the isothermal shock if there is one in the disk. On the other hand, the value 
of £ remains constant since the disk is inviscid. The energy transport rate per unit mass can 
be written as 

i e 1 2 1 2 GM 

e = - — + -vH a 2 , (1) 

2 r 2 2 7-1 r - r s ' K ' 

where v denotes the radial velocity (defined to be positive for inflow), r s = 2GM/c 2 , and the 
adiabatic sound speed is given by a = {^Pj p) 1 ! 2 for a gas with pressure P, mass density p, 
and specific heat ratio 7. Note that e is defined to be positive for an inward flow of energy 
into the black hole. Furthermore, in this semi-classical approach, e does not include the 
rest-mass contribution to the energy. 

Dissipation and radiative losses are unimportant in inviscid, advection-dominated disks, 
and therefore the flow is in general adiabatic and isentropic, except at the shock location. 
In this situation, the "entropy parameter," 

KEMWr 3 / 2 (r-r s ) , (2) 

is conserved throughout the flow except at the location of a shock (Becker & Le 2003). When 
a shock is present, we shall use the subscripts "-" and "+" to refer to quantities measured just 
upstream and just downstream from the shock, respectively. In the isothermal shock model, 
K and e have smaller values in the post-shock region (K + , e+) compared with the pre-shock 
region (K_, e_). The entropy and energy lost from the disk are carried off by the outflowing 
high-energy particles, and the total entropy of the combined (disk/outflow) system is higher 
when a shock is present, as expected according to the second law of thermodynamics. By 
combining equations (1) and (2) and differentiating with respect to radius, we can obtain a 
"wind equation" of the form (1/f) (dv/dr) = N/D, where 



WB -4 + g^ + , Ds *L-*. (3) 

T-d ( r _ r - s y ^ _|_ ;y r ( r _ rg j 7 + I 



Critical points occur where the numerator N and the denominator D vanish simultaneously 
(see Figure lb for a specific example). The associated critical conditions can be combined 
with the energy equation (1) to express the critical velocity v c , the critical sound speed 
a c , and the critical radius r c as functions of e and £. The value of K is then obtained by 
substituting r c , v c , and a c into equation (2). In general, one obtains four solutions for the 
critical radius, denoted by r cl , r c2 , r c3 , and r c4 . Only r cl and r c3 are physically acceptable 
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sonic points, although the types of accretion flows that can pass through them are different. 
Specifically, r c3 is an X-type critical point, and therefore a smooth, global, shock-free solution 
always exists in which the flow is transonic at r c3 and then remains supersonic all the way 
to the event horizon. On the other hand, r cl is an a-type critical point, and therefore any 
accretion flow originating at a large distance that passes through this point must display a 
shock transition below r c \ in order to cross the event horizon (Abramowicz & Chakrabarti 
1990). 

The radius of the isothermal shock, denoted by r*, is determined self-consistently 
along with the structure of the disk by satisfying the velocity and energy jump conditions 
(Chakrabarti 1989a) 

where M- = i>_/a_ is the upstream Mach number at the shock location and i?* is the shock 
compression ratio. Due to the escape of energy at the isothermal shock, the energy transport 
rate e drops from the upstream value e_ to the downstream value e+, and consequently 
Ae < 0. The flow must therefore pass through a new inner critical point located at f c3 < r*, 
which is computed using the downstream value e+. This structure represents the transonic 
solution of interest here. The kinetic power lost from the disk at the isothermal shock 
location is related to observable parameters by writing 

L jet = Mc 2 Ae , (5) 

where L )et and M are the observed kinetic power and the accretion rate for a specific source, 
respectively. Hence we conclude that Ae = Lj et /(Mc 2 ). For given observational values of 
M, M, and L- )Ct) only certain combinations of the specific angular momentum t and the 
upstream specific energy transport rate e_ will yield the correct value for Ae. Hence we 
can view e_ as the fundamental free parameter and determine t using the energy constraint 
Ae = Lj Ct /(Mc 2 ) once M, M, and Lj et are specified for a particular source. The two sources 
of interest here are M87 and Sgr A*, which correspond to models A and B, respectively. 



3. PARTICLE TRANSPORT EQUATION 

The model for the diffusive particle transport employed here is similar to the scenario 
for the acceleration of cosmic rays first suggested by Blandford & Ostriker (1978), although 
in our situation the shock maintains a fixed location in the disk and the seed particles for 
the acceleration process are provided by the accretion flow itself. The energy of the injected 
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Fig. 1. — Results obtained using the M87 values for M, M, and Lj et . In panel (a) we plot 
k,q (dashed line) and (solid line) as functions of the upstream energy transport rate e_ 
such that the energy conservation condition N esc E esc = L )et is satisfied. In panel (b) the flow 
velocity v (r) (solids lines) and sound speed a(r)[2/(7 + 1)] 1//2 (dashed lines) are plotted as 
functions of the radius for the model A parameters, which are indicated by the letter "A" 
in panel (a). In panel (b), the thin lines represent the shocked solutions and the heavy lines 
denote the shock- free solutions. See the discussion in the text. 



- 6- 



seed particles is given by Eq = 0.002 erg, which corresponds to an injected Lorentz factor 
T = E /(m p c 2 ) ~ 1.3, where m p is the proton mass. Starting with values for e_ and £ 
associated with a particular source, we can compute the injection rate for the seed particles, 
iVo, by setting N E = Mc 2 Ae = L Jct (see eq. [5]), which ensures that the energy injection 
rate for the relativistic seed particles is equal to the energy loss rate for the background gas at 
the isothermal shock location. The particle injection rate is therefore given by Nq = Lj Ct /Eo, 
which represents one of the fundamental self-consistency requirements for the model. 

The particle Green's function / G (e, r) describes the energy and spatial distribution of 
the relativistic particles injected with energy Eq at the shock location r*. The associated 
number and energy densities of the relativistic particles are given by n r = J °° AtiE 2 f G dE 
and U r = J °° 4tiE 3 f G dE, respectively. Including the escape of the relativistic particles at 
the shock location, the energy moments I n = J °° 4irE n f G dE of the Green's function satisfy 
the one-dimensional steady-state transport equation (Le & Becker 2004) 



4- (AnrHF n ) = AnrH 
dr 



;2-n)W4 + «^ (r _ r<) _ cW _ rJ 



3 dr 4irr*H* 



(6) 



where c is the speed of light, is the disk half-thickness at the shock location, and the flux 
F n is defined by 

_ (n + l)y dl n 

E N = — In-K—, (7) 

which is obtained by integrating over the vertical structure of the disk. The terms on the 
right-hand side of equation (6) represent first-order Fermi acceleration, the particle source, 
and the escape of particles from the disk at the shock location, respectively. The terms 
on the right-hand side of equation (7) describe the diffusion and advection of particles, 
respectively, and k denotes the spatial diffusion coefficient. The dimensionless parameter A 
is a constant that represents the microphysical processes controlling the vertical escape of 
energetic particles from the disk in the vicinity of the shock, as discussed below. 

To close the system, we must also specify the radial variation of the spatial diffusion 
coefficient k in equation (7). Since the radial dependence of k is not known precisely in this 
physical situation, we shall adopt the general form 

K(r) = K v(r)r a (j r -lj , (8) 

where kq is a dimensionless positive constant. Close to the event horizon, inward advection 
at the speed of light must dominate over outward diffusion. Conversely, in the outer region, 
we expect that diffusion will dominate over advection as r — > oo. A careful analysis of the 
differential equation for the relativistic particle number density n r obtained by setting n = 2 
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in equation (6) reveals that these conditions are both satisfied if a > 1, and in our work we 
shall set a = 2. 

The spatial diffusion coefficient k is related to the magnetic mean free path via the usual 
expression k = cA mag /3. Analysis of the three-dimensional random walk of the escaping 
particles then yields 

where = + k + )/2. We need to specify the value of k q in order to compute k_ and 
k + using equation (8). Since A is independent of the particle energy, it follows that all of 
the relativistic particles have the same escape probability per unit time in the vicinity of 
the shock, and therefore the mean energy of the escaping particles is given by E esc = U^/n*, 
where n* = n r (r*) and U* = U r (r*) denote the relativistic particle number and energy 
densities at the shock location. Note that E esc is proportional to E and independent of N . 
To obtain a self-consistent result, the product of the calculated escape energy E csc and the 
particle escape rate N esc at the shock location must be equal to the total energy lost from 
the background in the dynamical model, and therefore we set N esc E esc = Mc 2 Ae = Lj et . 
This condition is not automatically satisfied in general and therefore it constrains the free 
parameters e_ and k as explained below. 



4. APPLICATIONS TO M87 AND Sgr A* 

In our numerical examples, the disk structures for M87 and Sgr A*are determined based 
on observational estimates for the black hole mass M, the mass accretion rate M, and the 
jet kinetic power L- ]Ct . From observations of M87, we obtain M ~ 3 x 10 9 M (e.g., Ford 
et al. 1994), M ~ 1.3 x KT 1 M yr" 1 (Reynolds et al. 1996), and L jet ~ 10 43 " 44 erg s" 1 
(Reynolds et al. 1996; Bicknell & Begelman 1996; Owen, Eilek, & Kassim 2000). We shall 
adopt the average value L- ]Ct = 5.5 x 10 43 erg s -1 . For Sgr A*, observations indicate that 
M ~ 2.6 x 10 6 M (e.g., Schodel et al. 2002) and M ~ 8.8 x 10" 7 M yr" 1 (e.g., Yuan, 
Markoff, & Falcke 2002; Quataert 2003). The kinetic luminosity of the jet in Sgr A*is rather 
uncertain (see, e.g., Yuan 2000; Yuan et al. 2002). In our work, we will adopt the value 
quoted by Falcke & Biermann (1999), L jet ~ 5 x 10 38 erg s _1 . 

Our goal here is to compute self-consistently the mean energy of the escaping particles, 
E csc , while matching the observational values of M, M, and L jet for a specific source. We 
shall utilize natural gravitational units, with GM = c = 1 and r s = 2, except as noted. 
Following Narayan, Kato, & Honma (1997), we assume approximate equipartition between 
the gas and magnetic pressures, and therefore we set 7 = 1.5. Once the values of M, M, and 



- 8- 



Lj et have been specified for the source, we select a value for the fundamental free parameter 
e_ and then compute t based on the energy conservation condition Ae = Lj et /(Mc 2 ) (see 
§ 2). The velocity profile v(r) is determined by integrating numerically the wind equation 
(1/f) (dv/dr) = N/D, where N and D are given by equation (3). The associated solution 
for the isothermal sound speed a(r) can then be computed using equation (2) based on the 
fact that K is constant in the adiabatic flow (except at the shock if there is one). 

Once the velocity profile v(r) has been obtained, the transport equation (6) can be solved 
numerically to determine the number and energy density distributions for the relativistic 
particles in the disk, n r (r) = 12(f) and U r (r) = h{ r )- The energy conservation requirement 
N CSC E CSC = Ljet then allows us to solve for the corresponding value of the diffusion coefficient 
constant n . Based on the observational values for M, M, and L jet associated with M87, in 
Figure la we plot the results obtained for k as a function of e_. In general, one obtains two 
roots for n if e_ < e max , where e max = 0.001527 for the M87 parameters and e max = 0.001229 
for the Sgr A*parameters. No self-consistent solutions exist if e_ > e max , and a single root 
for is obtained if e_ = e max (see Fig. la). For illustrative purposes, we shall develop two 
detailed models by setting e_ = e max for M87 and Sgr A*. 

The parameter values for the two models are e_ = 0.001527, £ = 3.1340, k = 0.02044, 
A>o = 2.75 x 10 46 s~\ re, = 0.427877, A = 0.0124, n* = 1.19 x 10 44 cm" 3 , U* = 1.42 x 
10 42 ergcm- 3 , N esc = 4.61 x 10 45 s" 1 ,^ = 0.0119 erg, r* = 21.654, e+ = -0.005746, 
r cl = 98.524, r c3 = 5.379, f c3 = 5.659, A4_ = 1.125, R* = 1.897, and H, = 11.544 for 
model A (M87). For model B (Sgr A*), we have e_ = 0.001229, £ = 3.1524, k = 0.02819, 
A> = 2.51 x 10 41 s-\ k, = 0.321414, A = 0.0158, = 1.93 x 10 39 cm- 3 , U* = 2.10 x 
10 37 ergcm- 3 , A> csc = 4.56 x lO^s" 1 , E csc = 0.0109 erg, r* = 15.583, e+ = -0.008749, 
r cl = 131.874, r c3 = 5.329, f c3 = 5.723, A4_ = 1.146, R* = 1.970, and H* = 7.672. The 
corresponding dynamical solutions for v(r) and a(r) (with and without a shock) are plotted 
in Figure lb for model A. 

The terminal (asymptotic) Lorentz factor of the jet can be computed using = 
E csc /(m p c 2 ). In Figure la we plot as a function of e_ for the M87 parameters. Since 
in general we obtain two roots for k for each value of e_, we also obtain two values, 
with the larger one corresponding to the smaller kq root. For models A (M87) and B 
(Sgr A*), we obtain = 7.92 and = 7.26, respectively. The associated values for 
the mean energy boost ratio, E esc /Eo, are given by E esc /Eo = 5.95 and E esc /Eo = 5.45, 
respectively. Conversely, we find that in the shock- free models with the same values for e_, 
£, and k , the energy is boosted by only a factor of ~ 1.4 — 1.5. This clearly establishes 
the essential role of the shock in efficiently accelerating particles up to very high energies, 
well above the energy required to escape from the disk. The self-consistency conditions 
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NoE = Lj et and N esc E esc = L ]Ct can be combined to conclude that N csc = N Q E / E csc , which 
implies that < 20% of the injected particles escape (vertically) through the surface of the 
disk. The remainder of the particles are either advected into the black hole or else diffuse 
outward (horizontally) through the disk. Our predictions for the shock/jet location and 
the asymptotic Lorentz factor agree rather well with the observations of M87 reported by 
Biretta, Junor, & Livio (2002) and Biretta, Sparks, & Macchetto (1999). The shock location 
predicted by the model in the case of Sgr A*is fairly close to the value suggested by Yuan 
(2000). However, no reliable observational estimate for the jet Lorentz factor in Sgr A*is 
currently available. 

5. CONCLUSIONS 

In this Letter, we have demonstrated for the first time the physical consequences of a 
standing shock in an inviscid ADAF accretion disk for the acceleration of relativistic particles 
and the production of outflows. The global model presented here provides a single, coherent 
explanation for the disk/outflow structure based on the well-understood concept of first- 
order Fermi acceleration in shock waves. The theory is based on an exact mathematical 
approach to the solution of the combined hydrodynamical and particle transport equations. 
Our work may help to explain the observational fact that the brightest X-ray AGNs do not 
possess strong outflows, whereas the sources with low X-ray luminosities but high levels of 
radio emission do. We suggest that the gas in the luminous X-ray sources is too dense to 
allow efficient Fermi acceleration of a relativistic particle population, and therefore in these 
systems, the gas simply heats as it crosses the shock. Conversely, in the tenuous ADAF 
accretion flows studied here, the relativistic particles are able to avoid thermalization due to 
the long collisional mean free path, resulting in the development of a significant nonthermal 
component in the particle distribution which powers the jets and produces the strong radio 
emission. 

We have presented detailed results that confirm that the general properties of the jets 
observed in M87 and Sgr A*can be understood within the context of our disk/shock/outflow 
model. Although the specific numerical results discussed here are limited to the case of an 
inviscid disk, we believe that the inclusion of viscosity will not change the basic features 
of the inflow/outflow solutions obtained here because efficient first-order Fermi acceleration 
will also occur in viscous disks, provided they contain shocks. In addition to viscosity, future 
work must also address several other important considerations such as the effects of general 
relativity and radiative transport on the disk/shock structure. 

The authors are grateful to the anonymous referee for several useful suggestions that 
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